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Parametric optimization of neutral linear system 
with two delays with P-controller 


JOZEF DUDA 


In the paper the problem of parametric optimization of a linear neutral system with two 
delays accomplished with respect to the integral quadratic performance index is formulated and 
solved. A method of computing the value of performance index relies on determining such a 
Lyapunov functional defined on the state space that its value for the initial state is equal to that 
of performance index. In the paper the form of Lyapunov functional is assumed and a method 
of computing its coefficients is given. An example illustrating the application of the theory 
discussed is presented. It concerns a system with the P-controller designed to control a plant 
with two delays both retarded and neutral type. For such a system the value of the considered 
performance index is determined. 
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1. Introduction 


Lyapunov quadratic functionals are used to test the stability of systems, in compu- 
tation of the critical delay values for time delay systems, in computation of the expo- 
nential estimates for the solutions of time delay systems, in calculation of the robustness 
bounds for uncertain time delay systems, to calculation of a quadratic performance in- 
dex of quality for the process of parametric optimization for time delay systems. We 
construct the Lyapunov functionals for the system with time delay with a given time 
derivative. For the first time such Lyapunov functional was introduced by Repin [13] 
for the case of retarded time delay linear systems with one delay. Repin [13] delivered 
also the procedure for determination of coefficients of the functional. Duda [1] used the 
Lyapunov functional, which was proposed by Repin, for the calculation of the value of 
a quadratic performance index of quality in the process of parametric optimization for 
systems with time delay of retarded type and extended the results to the case of neutral 
type time delay system in [2]. Duda [3] presented a method of determining the Lyapunov 
functional for linear dynamic system with two lumped retarded type time delays in the 
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general case with no-commensurate delays and presented a special case with commen- 
surate delays in which the Lyapunov functional could be determined by solving of a set 
of ordinary differential equations. Duda [4] presented a method of determining the Lya- 
punov functional for linear dynamic system with two delays both retarded and neutral 
type time delay and in the paper [5] presented a method of determining the Lyapunov 
quadratic functional for linear time-invariant system with k-non- commensurate neutral 
type time delays. Infante and Castelan [7] construction of the Lyapunov functional is 
based on a solution of a matrix differential-difference equation on a finite time interval. 
This solution satisfies symmetry and boundary conditions. Kharitonov and Zhabko [11] 
extended the results of Infante and Castelan [7] and proposed a procedure of construc- 
tion of quadratic functionals for linear retarded type delay systems which could be used 
for the robust stability analysis of time delay systems. This functional was expressed by 
means of Lyapunov matrix, which depended on the fundamental matrix of time delay 
system. Kharitonov [8] extended some basic results obtained for the case of retarded 
type time delay systems to the case of neutral type time delay systems, and in [9] to 
the neutral type time delay systems with discrete and distributed delay. Kharitonov and 
Plischke [10] formulated the necessary and sufficient conditions for the existence and 
uniqueness of the delay Lyapunov matrix for the case of retarded system with one delay. 

The paper presents the parametric optimization problem for a system with two de- 
lays both retarded and neutral type time delay with P-controller. To the best of author’s 
knowledge, such parametric optimization problem has not been reported in the literature. 
An example illustrating the method is also presented. 


2. Formulation of the parametric optimization problem 


Let us consider the linear system with two delays both retarded and neutral type, 
whose dynamics is described by equations 


ax) pE = x(t) + Bx(t— 1) + Cu(t—r) 
0) = xo 


o +8) = (0) 
) = —Kx(t) 


x(t 
(1) 

x(t 

u(t 

t > to, 80 E |[-r,0), r>t>0 


A,B,DER™", x(t) ER", &eW'([-r,0),R") 


CER”®™,KEeR™”, ult) €R” 
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W'?((—7,0),R”) is a space of all absolutely continuous functions defined on the interval 
[—r,0) with values in R” which derivative is in L?(|—r,0), R”) -a space of a Lebesgue 
square integrable functions on interval [—r,0) with values in R”. 
We can reshape the equation (1) to the form 
a) pI — x(t) + Bx(t —t) —CKx(t —r) 
x(to) = Xo (2) 
x(to +8) = B(8) 


We introduce a new function y, defined by term 


y(t) =x(t)-—Dx(t—t) fort > to (3) 
Thus the equation (2) takes a form 


DO — ay(t) + (AD +B) x(t -1)—CKx(t — r) 
y(t) = x(t) — Dx(t = 1) 

y(to) = xo — D®(—T) 

x(to + 8) = (0) 


The state of the system (4) is a vector 


(4) 


S(t) = | yt) | for t >to (5) 


where x, €W!7?([—r,0),R")  x(0) = x(t +0) for @ € [—r,0) 


The state space is defined by the formula 


X =R" x W? ([-r,0), R”) (6) 


We search for such matrix K whose minimizes the integral quadratic performance 
index 


J= f y Eye)dt (7) 


3. The method of determining the value of the performance index 


On the state space X we define a Lyapunov functional, positively defined, differen- 
tiable, with the derivative computed on the trajectory of the system (4) being negatively 
defined. 
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V(S(t)) =y" (t)ay(t) + [ro \x(t-+0)d0-4+ (8) 
0 0 
+] | x"(t+0)8(0,c)x(t+0)dod0 
i 


for t > tọ where a = a? € R”, B e C!([-7,0],R"™”), 5 € C'(Q,R™"), Q = 
{(8,<) : 8 € [—r,0], ¢ € [0,0]}. C! is a space of continuous functions with continuous 
derivative. 

We identify the coefficients of the functional (8), assuming that its derivative com- 
puted on the trajectory of the system (4) satisfies the relationship 


ws) = —y" (t)y(t) fort >to (9) 


When the relationship (9) holds , we can easily determine the value of a square 
indicator of the quality of the parametric optimization, because 


J= fy" rar =V (S) (10) 


4. Determination the coefficients of the functional (8) 


We compute the derivative of the functional (8) on the trajectory of the system (4) 
according to the formula 


wis) Lea graa (s0) 39 


The derivative of the functional (8), calculated on the basis of the formula (11), is 
given by the formula 


dV (S(t)) 
dt 


fort to (11) 


T 
BOBO, 


=y" (t) lar + OA + 


+y" (r) Ra (B +AD) + B(0)D] x(t —1) +y" (t) [-20CK — B(—r)]x(¢-r)4 


0 


+ J yT (t) aeo- 4 ape tO) 30, o) x(t +0)d0+ 
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0 
$ f x" (t —1)[(B+AD)! B(0) + D8" (0,0)]x(t +0)d0+ 


0 
+ f x! (t —r) [-K™C"B(®) — 8(—r,6)] x(t +0)d0+ 


_f fate -o [EE (0,0) | age) erodi i 


-r 
From equations (12) and (9) we obtain the system of equations (13) to (19) 


aTa aa + BO +BY) -I (13) 

T 2 =s 
2a(B+AD)+B(0)D=0 (14) 
—2aCK — B(—r) =0 (15) 
A™B(@) — r OO) +30, 0) =0 (16) 
(B+AD)'B(8) + D7 8 (6,0) =0 (17) 
—K'C'B(0) —8(—r,0) =0 (18) 

05(0,6) , 05(6,6) _ 
Ə ` ð = a) 
for 8 € [—r,0], o € [—7,0]. 
From equation (14) it results that 

B(0) = —2a(B+AD) D7! (20) 


We put (20) into (13). After some calculations we get 


aP+P?o=I1 (21) 


where 


P = BD! (22) 
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From equation (21) we obtain the matrix a. 
From equation (17) we obtain 


ô” (0,0) = —A" B(0) — P”B(6) 
We put (23) into (16). We get 


aB(6) _ 
=P TBE) 


(23) 


(24) 


We get the initial condition for the differential equation (24) from the formula (15) 


B(—r) = —20CK 


(25) 


The solution of the differential equation (24) with the initial condition (25) is given 


by formula 


B(8) = —2exp ( P7 (04 r)) oCK 
for 6 € [—r,0]. 
The solution of the equation (19) is as below 
5(8,0) = 9(8—o) 
where ọ € C! ((—r,7],R"*"). 


From equations (27) and (18) we obtain 


5(—r,0) = ọ(—r — 80) = —KTCTR(0) 


for 6 € [-r,0]. 
Hence 
PÉ) = —K7C"B(—-§ —r) 
for é € [—7,0]. 


Taking into account (27) and (29), we get the formula 


5(0,0) = —K’C’B(o— 0—r) 


Upon taking the relation (26) into account, we get the formula 


5(8,0) = 2K'C' exp (—P" (o — 0)) aCK 


In this way we obtained all parameters of the Lyapunov functional (8). 


(26) 


(27) 


(28) 


(29) 


(30) 


(31) 
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5. Determination the value of the performance index 


According to the formula (10) the value of the performance index is given by the 


formula 
J =V (y(to),®) = y" (to) ay(to)+ (32) 
o 00 
+ i Y7 (to) B(0)®(0)d0-+ f f &! (8)3(0,0)(s)dode 
We put the sete (26) and (31) ae and we get 
(33) 


0 
J =y" (to) ay(to) —2 f Y” (to) exp (-P" (0+r)) aCK®(0)d0-+ 


0 0 
+2 / J DT (6)K™CT exp (PT (6 — 0)) 0CK®(o)dod6 
J 


6. Example 


Let us consider the system described by equation 


eo — ge = ax(t) + bx(t — T) — ckx(t — r) 

X(t) =x ER (34) 

x(to +0) = (6) 

t>t,x(t)ER, OE [-70), a,b,c,d,kER, r>t>0 
We can reshape the equation (34) to the form 
DU — ay(t) + (b+ ad) x(t — 1) — ckx(t — r) 
y(t) =x(0) -dalr =) es 
y(to) = x9 — d®(—T) 
x(to +8) = (0) 
a,b,c,d,kER, r>t>0 


Oe [-r,0), 


t > to, x(t) ER, 
We search for such parameter k whose minimizes the integral quadratic performance 


index 
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J= [7 OOd =V ((00).®) (36) 
The Lyapunov functional V is defined by the formula 


0 0 


V(S(t)) = ay?( + [von x(t +8) yao ff x(t +8) 8(0,0)x(t-+6)ded® (37) 


where according to (21), @ is given by formula 


1 d 
za (38) 
According to (26), the coefficient B is given by formula 
B(0) = —2ackexp ( ner 2) = a exp ( er 2) (39) 
According to (31), the element 6 is given by formula 
5(0,0) = 20c*k* exp ( Acia 2) = eae exp ( ee 2) (40) 
d b d 
The value of the performance index is given by formula 
0 
I= Ato) - Pro) | 2(8)exp (2S =") ae (ai) 


ce / j &(0)(c) exp (- n =) dod0 


The performance index is a quadratic function with respect to the variable k. 
To get the optimal value of the gain k we compute the derivative of the performance 
index with respect to k and we equal it to zero. 


dt (k) 
dk 

sy) _ 158008) 
2c  (0(8)@(c) exp CE dod® 


=0 (42) 


kopt = (43) 
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7. Conclusions 


The paper presents the parametric optimization problem for a system with two delays 
both retarded and neutral type time delay with P-controller.A method of computing the 
value of performance index relies on determining such a Lyapunov functional defined 
on the state space that its value for the initial state is equal to that of performance index. 
In the paper the form of Lyapunov functional is assumed and a method of computing its 
coefficients is given. 
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